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ABSTRACT 

VV124 (UGC 4879) is an isolated, dwarf irregular/dwarf spheroidal (dlrr/dSph) transition-type 
galaxy at a distance of 1.36 Mpc. Previous low- resolution spectroscopy yielded inconsistent radial 
velocities for different components of the galaxy, and photometry hinted at the presence of a stellar 
disk. In order to quantify the stellar dynamics, we observed individual red giants in VV124 with the 
Keck/DEIMOS spectrograph. We validated members based on their positions in the color-magnitude 
diagram, radial velocities, and spectral features. Our sample contains 67 members. The average 
radial velocity is (v r ) — —29.1 ± 1.3 km s , in agreement with the previous radio measurements 
of H I gas. The velocity distribution is Gaussian, indicating that VV124 is supported primarily by 
velocity dispersion inside a radius of 1.5 kpc. Outside that radius, our measurements provide only an 
upper limit of 8.6 km s" 1 on any rotation in the photometric disk-like feature. The velocity dispersion 
is g v = 9.4 ± 1.0 km s -1 , from which we inferred a mass of Mi/ 2 = (2.1 ± 0.2) x 10 7 Mq and a mass- 
to-light ratio of (M/Lv)i/2 = 5.2 ± 1.1 Mq/Lq, both measured within the half-light radius. Thus, 
VV124 contains dark matter. We also measured the metallicity distribution from neutral iron lines. 
The average metallicity, ([Fe/H]) = —1.14 ± 0.06, is consistent with the mass-metallicity relation 
defined by dwarf spheroidal galaxies. The dynamics and metallicity distribution of VV124 appear 
similar to dSphs of similar stellar mass. 

Subject headings: galaxies: individual (VV124) — galaxies: dwarf — Local Group — galaxies: kine- 
matics and dynamics — galaxies: abundances 



1. INTRODUCTION 

The dwarf galaxy VV124 is one of the "newest" mem- 
bers of the local univers e. The discovery of VV124 
was first published in the IVorontsov-Velvaminovl ([19590 
Atlas and Catalog of Interacting Galaxies, though the 
galaxy does not have an interacting neighbor. The 
CfA galaxy redshift survey (Huchra ct al. 1983) reported 
its redshift as cz = 600 ± 100 km s _1 . In the Hub- 
ble flow, this redshift corresponds to a distance of ~ 
10 Mpc. At that distance, VV 124 would have be- 
longed to a group of th ree galaxies ([Turner fc Gottlll976l : 
iGeller fc Huchralll983l) . However. IKopylov et ^ (2008) 
noticed resolved stars in Sloan D ig ital Sky Survey images 
(jAdelman-McCarthv et al.ll2007h . IKopylov et alJ imaged 
VV124 to a depth of V — 25.6, and they confirmed that 
it is in fact much nearer than 10 Mpc. They measured a 
distance of 1.1 ±0.1 Mpc based on the tip of the red giant 
branch (TRGB). This new distance placed the galaxy on 
the periphery of the Local Group. Despite its proxim- 
ity to the group, VV124 is extremely isolated. Its nearest 
neigh bor, Leo A, is anoth er isolated dwarf galaxy 700 kpc 
away ([Jacobs et al.ll201ll ). 

The revision to VV124's distance drastically affected 
its inferred star formation history. From low-resolution 
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spectroscopy, Uames et alJ ([2004) measured its Ha star 
formation rate as (5.0 ± 2.4) x 10~ 3 M yr _1 , but that 
rate was based on a distance of 1 0.5 Mpc. Correcte d 
to the true distance of 1.36 Mpc ([Jacobs et al.l 12011). 
the rate dropped to 8.4 x 10~ 5 M w yr _1 (a point also 
noted by iBellazzini et all 1201 lal hereafter Blla). The 
new distance also inspire d deep photometry in addi- 
tion to that obtained bv IKopylov et alJ (|2008[ ). iBllal 
obtained wide-field g, r photometry of VV124 at the 
Large Binocular Telescope (L BT) to study the s truc- 
ture of the galaxy in detail. Uacobs et alJ ([201 1|) ob- 
tained even deeper V and I photometry over a smaller 
field with the Hubble Space Telescope (HST) Advanced 
Camera for Surveys. They determined that 93% of the 
stars a re older than 10 Gyr. T he rest are younger than 
1 Gyr. IBellazzini et alJ ([2011bfi used the same HST im- 
ages as I Jacobs et al.l to show that the oldest, most metal- 
poor population is significantly more extended than the 
younger st ars. Th e galaxy-averaged star formation rates 
that Uacobs et al.1 derived were (3.0±0.2) x 10 -4 M© yr" 1 
over the last 0.5 Gyr and (8.0 ± 0.5) x 10~ 4 M© yr" 1 
over the last 1 Gyr. Together, the HST photometry and 
Ha spectroscopy indicated that VV124 is a typically old 
dwarf galaxy that experienced a small burst of star for- 

mati on ~ 1 Gyr a go. 

Tik honov et al.l (|2010D expanded IKopylov et all 's 
(2008) study. They classified VV124 as transition- 
type galaxy, caught between being an dlrr and a dwarf 
spheroidal galaxy (dSph). They noticed that the young, 
blue stars lie in a thin, disk-like structure, whereas the 
older, red giants occupy a thicker configuration. How- 
ever, the small number of very young stars was incon- 
sistent with an irregular-only classification. Therefore, 
they concluded that VV124 is a rare dlrr/dSph galaxy. 
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iGrebel et al.l (J2003D listed only five other such mem- 
bers of the Local Group: Antlia, Aquarius, Phoenix, 
Pisces, and KKR 25. This class of galaxies is espe- 
cially interesting because their existence supports the 
theory that dSphs are dlrrs that lost their gas through 
ram pres sure stripping and ti da l interaction s with other 
galaxies (jMaver et al.l 12001 alibi . 120061 |200l . There is 
strong observational support for the removal of gas from 
dwar f galaxies in dense environm ents (e.g., iGeha et ah! 
120061: iGrcevich fc Putmanl I2009T ) . and the similarity of 
the star formation histories except in the past Gyr of 
dlrrs, dSphs, and dlrr/dSphs fu rther supports the trans- 
formation of dlrrs into dSphs ()Weisz et al.l 1201 ID . Al- 
though there are counter examples, such as the compara- 
tively young Leo A dlrr (|Cole et al.l 120071 ). the star for- 
mation history of VV124 less recent than 1 Gyr appears 
much like other dlrrs and dSphs. Except for the recent 
burst. iJacobs et all (|2011| ) reported that almost all of its 
stars are ancient, like nearly all dSphs. 

The dynamics of VV124 have been unclear. Despite 
the corrected redshift, measurements o f the galaxy's 
radial velocit y have been inc on sisten t. iKopvlov et al.l 
(2008) and iTikhonov et~aT1 (|2010h obtained low- 
resolution spectroscopy of individual stars, resolved H II 
regions, diffuse stellar light, and diffuse emission. The ve- 
locities varied widely, but the y seemed to cluster around 
v r = —70 km s _1 . Blla measured the H I veloc- 
ity structure with the Westerbork Synthesis Radio Tele- 
scope. They measured a heliocentric radial velocity of 
v r = —25 ± 4 km s _1 with a dispersion of ~ 11 km s _1 
and a full range of about —45 to km s _1 . iBllal also 
measured the velocities of a blue supergiant and a young 
star cluster, but the cluster velocity was inconsistent with 
(more negative than) the H I velocity by 2.7<r. So far, 
the discrepancy between the velocities of the neutral gas 
a nd you nger populations has not been resolved. 

IBllal also observed what appears to be a stellar disk. 
VV124 is elongated (e = 0.44±0.04, PA = 84°±11°), but 
its shape cannot be described by a single ellipsoid. There 
is an even more elongated component whose major axis 
is aligned with th e main spheroid. The second compo- 
nent, which lBl lal called "wings," resembles a stellar disk, 
which should r otate coherently, unlike the spheroid. As 
an alternative, IBllal also proposed that the wings might 
be tidal tails instead of a disk. However, VV124's ex- 
treme isolation makes strong tidal interaction unlikely. 

In order to provide closure to the discrepancies in ve- 
locity measurements and to test for rotation, we observed 
red giants in VV124 with the Deep Extraga lactic Imag- 
ing M ulti-Object Spectrograph (DEIMOS, iFaber et all 
2003) on the Keck II telescope. The median veloc- 
ity uncertainty of the individual stellar measurements is 
4.1 km s , which is small enough to test for the presence 
of dark matter. The spectra are also of sufficient qual- 
ity to measure metallicities from neutral iron lines. The 
median uncertainty in metallicity is 0.21 dex. This arti- 
cle presents the results of our dynamics and metallicity 
studies. 

2. OBSERVATIONS AND REDUCTIONS 

2.1. Target selection 

Targets were selected from lBl nil 's photometric catalog. 
DEIMOS slitlets cannot overlap in the dispersion direc- 



tion, and the VV124 field is dense enough that we needed 
to prioritize the targets for placement on the DEIMOS 
slitmasks. The targets were prioritized by magnitude and 
physical location with respect to the center of VV124. 
Brighter targets were preferred for higher signal-to-noise 
ratios (SNRs). We defined X as the distance along the 
major axis and Y as distance along the minor axis. This 
is the same nomenclature introduced by IBllal 

We used the dsimulatoiQ program to select targets for 
placement on DEIMOS slitmasks. The program maxi- 
mizes the number of high-priority targets within user- 
defined constraints. We required slitlets to be at least 
4" long and separated by 0.35" for adequate sky sub- 
traction. Slit widths were 0.7", a choice that balanced 
observing efficiency with spectral resolution, which corre- 
sponds to precision in the radial velocity and metallicity 
measurements . 

We prioritized targets for spectroscopy in the following 
order: 

1. Brighter red giant branch (RGB) candidates near 
the major axis (V < 23.7 and \Y\ < 2') 

2. Fainter RGB candidates near the major axis (V > 
23.7 and \Y\ < 2') 

3. Bright, blue non-RGB stars near the major axis 
(V < 23.7, g - r < 1.1, and \Y\ < 2') 

4. RGB candidates away from the major axis (\Y\ > 
2') 

5. Any other non-RGB candidates 

RGB candidates are defined as r > 22 and 0.7 < g — r < 
1.4. 

One pair of slitmasks comprised our observations. The 
slitmasks, subtending roughly 16' x 4', were roughly cen- 
tered on the galaxy at position angles of 86° and —94° 
(i.e., rotated 180°). The target lists for the two slitmasks 
were almost identical. This ensured that we had two in- 
dependent measurements of the radial velocities for al- 
most all of the stars. The target list included 131 unique 
targets, of which 106 were included on both slitmasks. 

Figure [T] shows the color-magnitude diagram (CMD) 
for VV124. Spectroscopic targets are identified by col- 
ored points. Filled points indicate spectroscopically con- 
firmed member stars (see Section d]). The colors corre- 
spond to the priorities in the list above. Figure [2] shows 
the coordinates of the stars with the same color coding. 

2.2. Observations 

We observed the two slitmasks, called vvl24a and 
vvl24b, on 2011 January 30 (UT). Table Q] lists the ex- 
posure times, seeing, number of targets, and number of 
VV124 members for each slitmask. The following night 
was also allocated to this project, but high humidity and 
fog on that night prevented any observing of VV124. 

We used the 1200 lines mm -1 grating centered at 
7800 A. This grating, combined with the slit widths of 
0.7" , yielded a line spread profile with a FWHM of 1.2 A. 
The resolving power is about R ~ 7000 near the Ca II 

5 http : //www . ucolick . org/$\sim$phillips/deimos_ref /masks . html 
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Table 1 

DEIMOS Observations 



Mask Seeing (") Targets Members Total Exp. Time (s) Individual Exp. Times (s) 

vvl24a 0.57-0.87 121 61 13200 

vvl24b 0.55-0.65 120 64 13600 



6 x 1800, 2 x 1200 
6 x 1800, 1600, 1200 
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Figure 1. CMD from LBT photometry (Blla). Colored points 
identify stars selected for DEIMOS spectroscopy. The color signi- 
fies the priority of the star for spectroscopic selection (Section 12.10 . 
Filled points are spectroscopically confirmed member stars, and 
hollow points are non-members. The right axis gives the approxi- 
mate Cousins V magnitude. 
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Figure 2. DEIMOS slitmask footprints overlaid on the LBT as- 
trometric catalog (Blla). The colors and types of points have the 
same meanings as in Figure[TJ The axis labels give the displacement 

from the center of VV124 (a = 9 h 16 m 03 s , <5 = +52°50'31"). The 
top and right axes give this displacement in kpc for an assumed dis- 
tance of 1.36 Mpc (Jacobs ct al. 2011). The positions and shapes 
of the two DEIMOS slitmasks are outlined. 

triplet (CaT) at 8500 A. We obtained a 1 sec exposure of 
internal Ne, Ar, Kr, and Xe arc lamps for wavelength cal- 
ibration and three 12 sec exposures of an internal quartz 
lamp for flat fielding. Although DEIMOS is mounted on 
the Nasmyth platform of Keck II, it rotates to track the 
sky. An active flexure compensation system maintains 
stability to within 0.2 pixels, corresponding to 3 /xm on 
the focal plane, which is equivalent to 0.02" in the spatial 
direction and to 0.07 A in the dispersion direction. 
Individual exposure times for the VV124 slitmasks 



ranged from 20 to 30 minutes. Most of the exposures 
were 30 minutes. We checked the alignment of the slit- 
mask about once per hour by placing the grating in ze- 
roth order. In this direct imaging mode, bright alignment 
stars aligned with 4" boxes in the slitmask. We mea- 
sured the seeing by fitting Gaussian profiles to the point 
spread functions of the alignment stars. Seeing remained 
low and stable for the entire night (see Table [I]). 

2.3. Reductions 

We reduced the raw images into one-dimensional 
spectra with the spec2c0 software developed by the 
Deep Extragalact ic Evolutionary Probe 2 team (DEEP2, 
iDavis et al.ll2003h . The software excised the footprint of 
each spectrally dispersed slitlet from each exposure of 
the flat field, arc lamps, and science targets. Because 
the DEIMOS focal plane is curved, the program recti- 
fied the slitlet images by tracing and straightening both 
the edges of the flat field and the roughly orthogonal 
arc lines. Then, the program generated a flat field im- 
age from the three exposures of the quartz lamp. This 
flat field was used to remove the slit response function, 
dust, and cosmetic imperfections in the science images. 
The spec2d software also fitted a fifth-order Legendre 
polynomial in wavelength to the pixel positions of arc 
lines. For each exposure of each slitlet, spec2d com- 
puted the sky emission as a function of wavelength and 
subtracted it from the two-dimensional spectrum. All of 
the science exposures for each slitmask were averaged 
with inverse variance weighting and cosmic ray rejec- 
tion. Finally, the object spect rum was extracted with 
optimal weighting (Home 1986). The final result was a 
clean, one-dimensional spectrum, wavelength-calibrated 
to a precision of 0.01 A. The spectrum was not flux 
calibrated, and it did not need to be for our purposes. 

Some spectra suffered from artifacts or such low SNR 
that measurement of the radial velocity was not possible. 
There were 13 such spectra. They are included and iden- 
tified in Table [21 but we do not consider them further. 

The spec2d software also computed the error spectrum 
from Poisson statistics of the raw frame. The errors were 
adjusted appropriately during flat fielding, and the one- 
dimensional error spectrum was extracted identically to 
the flux spectrum. We used the error spectra for estima- 
tion of uncertainties in radial velocities (Section [3]) and 
metallicities (Section [6]). 

As mentioned in Section 12. 1\ the two slitmasks in- 
cluded 106 duplicate targets. We measured the precision 
of our velocity measurements using the independently 
measured radial velocities from the two slitmasks. How- 
ever, for the final list of radial velocities, we used a com- 
bined spectrum. Wc summed the two onc-dimcnsional 
spectra for each duplicate object with inverse variance 

6 http : //deep . berkeley . edu/spec2d/ 
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weighting. 

We removed telluric absorption using archived 
DEIMOS spectra of hot stars. We used the s ame hot star 
spectr a and spectral division procedure as iKirbv et al.1 
(2008). However, our science exposures took place over 
an entire night. The varying airmass and water vapor 
made it impossible to remove the telluric absorption ac- 
curately. Regions of moderately strong telluric absorp- 
tion were not considered for metallicity measurements. 

The metallicity measurements discussed in Section |6] 
required normalization of the spectral continuum. We 
followed t he continuum norm alization procedure de- 
scribed by IKirbv et all (|2009t their Section 3.4). Af- 
ter normalization, we calculated the SNR of each spec- 
trum. First, we computed the median absolute devia- 
tion (m.a.d.) from the co ntinuum of pixels in "contin- 
uum regions" (defined by IKirbv et al.l 120081 ) in the rest 
wavelength ranges 7400-7500 A, 7750-8100 A, and 8400- 
8900 A. These selections excluded most telluric absorp- 
tion and TiO bands. Then, we removed all pixels that ex- 
ceeded five times the m.a.d. and re-calculated the m.a.d. 
with this slightly trimmed set of pixels. The result was 
the SNR per pixel. To convert to SNR per A, we mul- 
tiplied by ^/l/0.33, where 0.33 is the pixel scale in A 
per pixel. This measurement of SNR is sensitive to in- 
accuracies in the continuum placement, which will ar- 
tificially decrease the SNR. However, errors in contin- 
uum placement dominate the SNR calculation only at 

SNR > 60 A , which is much higher than any member 
star in our sample. 

3. RADIAL VELOCITY MEASUREMENTS 

We measured radial velocitie s and their un c ertain - 
ties following the procedure of iSimon fe Gehal (|2007| ). 
First, we measured the radial velocity of the star by 
cross-correlating the observed spectrum with a library 
of template stars and galaxies observed with DEIMOS. 
The template s pectra were the same as those used by 
ISimon fc Gehal They were observed with a nearly iden- 
tical configuration of DEIMOS and therefore have the 
same resolution as our science spectra. The strong atmo- 
spheric A and B bands (7588-7700 A and 6862-6950 A) 
were excluded from the cross-correlation. Three weaker 
telluric bands (7167-7315 A, 8210-8325 A, and 8900- 
9200 A) were given weights of 10% of other pixels in the 
cross-correlation. We adopted the velocity correspond- 
ing to the cross-correlation peak of the template with 
the lowest x 2 when compared to the observed spectrum. 
We call this velocity u bj ■ Figure [3] shows typical spec- 
tra of VV124 member stars along with the best- fitting 
template spectra. 

Next, we computed a correction for mis-centering of 
the star in the slitlet ()Sohn et al.1 120071 ). For this step, 
we used the observed spectrum before correction for tel- 
luric absorption. We cross-correlated this spectrum with 
the spectrum of HR 1641, a hot, rapidly rotating star. 
Only regions subject to signif icant telluric absorption 
(the same regions used by ISimon fe Gehal [2007) were 
considered in the cross-correlation. We call the veloc- 
ity obtained from this cross-correlation Utell- The radial 
velocity of the star is v r = w bj — v tc \i — fhciio, where 
^hoiio is the heliocentric velocity correction appropriate 
for Keck Observatory on the observation date. 
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Figure 3. DEIMOS spectra of three red giant members of VV124. 
Only a small portion of the spectrum around the CaT is shown. 
For display, the spectra have been normalized to the continuum 
and smoothed with a Gaussian kernel (a = 0,7 A). The rad ial 
velocity template spectrum, provided by Simon & Gcha (2007), is 
shown in red in all three panels. The top and bottom panels show 
the spectra with the highest and lowest SNRs in our sample. The 
middle panel shows the faintest member star in our sample. 



We estimated uncertainties in v r by Monte Carlo re- 
sampling of the spectra. For each of 1000 realizations, we 
added Gaussian random noise to the spectrum. The dis- 
persion of the Gaussian probability distribution at each 
pixel was equal to the estimated Poisson error of that 
pixel. Then, we recomputed u bj and t>t e ii for each real- 
ization. For computational efficiency, we considered only 
the template spectrum that best matched the original 
spectrum. The standard deviation of the 1000 measure- 
ments of v r is itmc • 

From duplica te v r measurements of the same stars, 
ISimon fc Gehal realized that ctmc was an incomplete es- 
timate of the radial velocity uncertainty. In order to 
explain the dispersion among the repeat measurements, 
they added a systematic error term, <7 sys , in quadrature 
with <tmc- The systematic term was chosen such that 
the radial velocity differences divided by their total un- 
certainties were distributed normally with unit standard 
deviation. 
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Because we observed 106 duplicate objects on both 
DEIMOS slitmasks, we also calculated a systematic er- 
ror term. We excluded galaxies, obvious foreground 
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Figure 4. Radial velocity measurement differences of the same 
stars observed on two different slitmasks. The differences are di- 
vided by the quadrature sum of the measurement uncertainties, 
which include noise and systematic error terms. The solid curve 
is a Gaussian with unit standard deviation. The systematic error 
term, <7 S ys, was adjusted to ensure that the standard deviation is 
unity. 
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Figure 5. Top: Absolute radial velocity differences of repeat mea- 
surements of the same stars on two different slitmasks versus SNR. 
These measurements were necessarily based on spectra with about 
half the exposure time of the final, stacked spectra. Bottom: Esti- 
mated error on v r as a function of SNR for all member stars. The 
total error is the quadrature sum of the random (Monte Carlo) and 
systematic error terms. 



dwarfs with Na I 8190 absorption, and other objects fail- 
ing any membership criterion other than radial velocity 
(see Section [4j. From the 57 remaining stars, we com- 



puted <7 SVS = 2. 21 km s^ 1 — which is the same value that 
iSimon feG cha derived — by solving Equation [1] Figure [4] 
shows the distribution of differences in repeat v r mea- 
surements divided by their estimated uncertainties. The 
width of the distribution matches a unit Gaussian be- 
cause er sys was tuned to force such a match. 

Figure [5] shows how velocity errors and uncertainties 
trend with SNR. The v r differences between repeat mea- 
surements (top panel) are based on exposures from one 
slitmask only. Therefore, their SNRs are about l/v2 
of the SNR of the stacked spectra used for remainder 
of our analysis. The estimated uncertainties (bottom 
panel) decrease with increasing SNR, as expected. Only 
the random error term, ctmCj decreases with SNR. The 
systematic error term, a sys , is constant. 

4. MEMBERSHIP 

We used several criteria for determining membership. 
First, we eliminated stars based on their positions in the 
CMD. Then, we eliminated stars with spectral features 
not found in red giants. Finally, we eliminated stars 
based on their radial velocities. 

4.1. Color-Magnitude Diagram 

Blla identified the extinction-corrected magnitude of 
the TRGB as r = 22.61 ± 0.06. In selecting spectro- 
scopic targets, we prioritized stars with r > 22, which 
conservatively allowed for large error in the TRGB mag- 
nitude. It also allowed some asymptotic giant branch 
(AGB) stars. In order to fill the slitmasks, we also in- 
cluded some targets with r < 22. 

Stars brighter than the TRGB are not likely members. 
We carefully examined all stars with tq < 22.49 (2a 
brighter than the TRGB). After excluding stars failing to 
meet the membership criteria in the following sections, 
four stars remained. Star 8116 showed broad CaT lines. 
The width and shape suggested pressure broadening typ- 
ical of dwarf stars. Therefore, we ruled 8116 as a non- 
member. The spectrum of star 8445 looked in all respects 
like a red giant. Although the star is slightly brighter and 
bluer than the TRGB, we kept 8445 in the member list 
under the premise that it is an AGB star. The spec- 
trum of star 11163 showed very strong cyanogen bands. 
If this carbon star is an AGB member of VV124, it could 
be significantly brighter and redder than the TRGB. In 
fact, that is its location in the CMD. For this reason, we 
kept 11163 in the membership list. Finally, star 13649 is 
solidly in the "wall" of foreground dwarfs in the CMD, 
and its Na I 8190 equivalent width (EW) is very close 
to the cut-off for membership (see the next section) . We 
ruled 13649 a foreground dwarf for the combination of its 
position in the CMD and moderately strong Na I 8190 
doublet. 

4.2. Spectral Features 

Inspection of the DEIMOS spectra revealed that many 
of our targets were background galaxies or foreground 
dwarf stars. First, we excluded all objects with emission 
lines or broad, redshifted Ca H and Ca K in absorption. 
These 20 objects are background galaxies. 

Second, we excluded stars with very strong, broad 
Na I 8190 doublets. The strengths of these lines are espe- 
cially sensitive to surface gravity. That sensitivity makes 
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Figure 6. Summed EW of the Na I 8190 doublet as a function 
of surface gravity and [Na/H], The summed EW exceeds 1 A only 
for dwarf stars (log g > 4.5). 



them excellent discriminators between giants and fore- 
ground dwarf s tars (Spinrad & Taylor 19 71[lCohen 1978; 
iSchiavon et~aT]|1997UGilbert et al 1 12006). We measured 
the EWs of each of the two lines in the doublet by fitting 
two independent Gaussians. We computed Monte Carlo 
uncertainties by resampling the spectrum 100 times. For 
each pixel in each resample, we added noise from a Gaus- 
sian random distribution with a standard deviation equal 
to the estimated Poisson noise of that pixel. We fitted 
Gaussians to each resample, and we adopted the stan- 
dard deviation among their EWs as the error on the EW. 
Table [5] gives the summed EWs and their uncertainties 
for stars where the doublet was measurable above the 
noise. 

We computed a threshold value for the summed EW 
of the doublet by looking up EWs in a grid of non- 
local t hermodynamic e quilibrium (NLTE) Na abun- 
dances (jLind et al.ll2011h . We sampled EWs from log g = 
0.5 to 4.7 and [Na/H] = 0.0 to -2.0. We found the ap- 
propriate effective temperature for each combination of 
surface gravity and abundance (assumi ng [Na/Fel = 0) 
from a 10 Gyr Yonsei-Yale isochrone (|Demarque et al.l 
120041) . Interpolations in the NLTE grid yielded Figure El 
Some points are missing from the figure because they are 
also missing from the NLTE grid. Regardless, it is clear 
that the doublet increases sharply in EW for dwarf stars 
with logg > 4. The summed EW of the doublet exceeds 
1 A for stars with [Na/H] < only at surface gravi- 
ties logg > 4.5. We identified 17 stars with a summed 
Na I 8190 EW greater than 1 A. These 17 stars are 
foreground dwarfs. 

Counting star 8116, excluded for its broad calcium 
lines and for being brighter than the TRGB, we elim- 
inated 38 objects based on spectral features. Another 
three stars showed strong CN bands. One of these, 
11163, is discussed in Section |4~T1 The other two stars, 
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Figure 7. Heliocentric radial velocity distribution. The solid 
curve is a Gaussian distribution corresponding to the measured 
mean velocity and velocity dispersion. The dashed lines enclose 
99% (±2.58<t) of the member stars. The figure label gives the av- 
erage heliocentric velocity and the velocity dispersion corrected for 
measurement uncertainty (see Section . Only stars passing all 
membership cuts except velocity are shown. 

9378 and 13043, pass all membership criteria, and we 
retained them as members. 

4.3. Radial Velocities 

The median radial velocity uncertainty in our sample 
of member stars is 4.1 km s" 1 , on the order of the veloc- 
ity dispersion of a small dSph. Therefore, it is important 
to consider velocity errors carefully when computing the 
velocity dispersion. [Walker et~a l. (2006) devised a for- 
malism to take these errors into account. The procedure 
involved maximizing the likelihood that the average ra- 
dial velocity, (v r ), and the intrinsic velocity dispersion, 
a v , accurately describe the observed velocity distribu- 
tion of candidate member stars in the presence of ob- 
servational error. The process was iterative. We chose 
(v r ) = —30 km s _1 and a v = 9 km s^ 1 as rough guesses 
for the first iteration. In each subsequent iteration, only 
stars within 2.58cr„ (99% of a normal distribution) of (v r ) 
were considered. In fact, just one iteration was required 
to converge on the final member sample and the final 
values of (v r ) and o~ v . The uncertainties on (v r ) and o~ v 
were determined from the likelihood covariance matrix. 
Table[3]gives the fiducial values for (v r ) and a v and their 
uncertainties. 

We used (v r ) and a v to determine the final membership 
cut. If a normal distribution is a good description of the 
velocity distribution, then the velocity range (u r )±2.58<7„ 
includes 99% of the members. This cut excluded 6 stars 
and retained 67 members. Some galaxy dynamics studies 
suggest that a cut of 3ct is a better choice. Adopting a 3c 
cut for VV124 would not have affected the membership 
of any star. 

Table [2] lists all of the spectroscopic targets. The table 
gives astrometric and photometric data, as well as quan- 
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tities derived from the spectra. The last two columns of 
the table indicate whether the object is a member star 
and the reasons that stars were considered non-members. 
Figure [7] presents the velocity distribution of VV124. 
The velocities are distributed normally. The skew- 
ness, a measure of asymmetry, of the v r distribution is 
0.14±0.29. The kurtosis, a measure of the deviation from 
Gaussianity, is — 0.26±0.58. Both parameters are consis- 
tent with zero, indicating conformance to a normal dis- 
tribution. The vertical, dashed lines show (v r ) ± 2.58(7,,, 
and the solid curve shows a perfect normal distribution 
(not a fit) for those parameters, normalized to the num- 
ber of member stars. 

5. DYNAMICS 

5.1. Mean Radial Velocity and Velocity Dispersion 

The mean radial velocity, (v r ) = —29.1 ± 1.3 km s _1 , 
agrees with the r adial v elocity of the H I gas, (v r ) = 
—25 ± 4 km s _1 (|Blla[ ). The agreement is especially 
reassuring because previous optical measurements of ra- 
dial velocities showe d signific antly more negativ e v r . 
iKopylov et all ((20081 ) and iTikhonov etldl (|2010l ) ob- 
served two blue supergiants: bll with v r ■ —90 i 



15 km s and bl2 with v r 



-82 ± 15 km s" 



They 



also estimated the average radial velocity of unresolved 
stellar light to be (v r ) = —70 ± 15 km s^ 1 . Radial 
velocities of resolved and diffuse H II gas, measured 
from the [O III] 50 07 emission line, ranged from —36 to 
—71 km s _1 . IBllal revised the velocity of the supergiant 
bl2 to v r = —44 ± 18 km s _1 . They also observed a star 
cluster, CI, with v r = —86 ± 20 km s _1 and diffuse Ha 
emission with (v r ) = —6 ± 30 km s . The individual 
sources bll, CI, and one of the H II regions have radial 
velocities more than la inconsistent with our velocity 
cuts for membership . The stellar i nteg rated light mea- 
surement of (v r ) by IKopylov et al.l and ITikhonov et al.l 
is 2.7a discrepant from our measurement. In every case, 
the discrepancies are such that the other measurements 
are more negative than ours. We are confident that our 
measurement of (v r ) is accurate because it agrees with 
the H I measurement. Moreover, our estimate is based 
on h igher-resolution spect r a and a much larger s ample 
than IKopylov et al.l (J2008D . ITikhonov et~aH (J2010D . and 
iBlTat 

There are two possible causes of the discrepancy. First, 
there could have been some systematic error in the 
measurements of v T from the low- resolution spectra of 
IKopylov et al.l a nd ITikhonov et al.l The upward revision 
of v r for bl2 by iBlla may support this conclusion. On 
the other hand, IBllal also measured a very negative v r 
for the star cluster CI. We also note that we excluded 
three stars with —75 km s _1 < v r < —55 km s _1 as 
non-members because their velocities are too negative. 
These velocities are similar to those of bll and CI. It 
may also be relevant that the H I gas southeast of the 
center of the galaxy h as mor e negative v r than the gas 
in the central regions (JBllal) . Thus, VV124 may have 
a non-equilibrium dynamical component consisting of a 
young stars and star-forming regions. This component 
would have more negative v r than the old red giants and 
neutral gas. The data available do not point to a defini- 
tive conclusion. Higher resolution observations of bll, 
CI, other young stars, and star-forming regions will bet- 



ter define VV124's velocity structure. 

Unre solved binari e s coul d inflate the velocity disper- 
sion. iMinor et al.l (|2010l ) concluded that this infla- 
tion is unlikely to be more than 30% in galaxies with 
4 kms" 1 < a v < 10 km s" 1 . In a separate study, 
iMcConnachie fc Cotel (|2010l ) arrived at a similar conclu- 
sion, but cautioned that even a modest binary fraction 
could cause an ultra-faint dSph-sized stellar population 
free of dark matter to appear to have the velocity disper- 
sion of a dark matter-dominated system. In a confirma- 
tion of the contrib ution of unreso l ved b inaries to mass 
estimates of dSphs. lKoposov et al.l (|2011|) revised the ve- 
locity dispersion of the Bootes I dSph f rom > 6 km s _1 
(|Munoz et al.l 120061: iMartin et ail 120071 ) to two compo- 
nents: one population at 2.4 km s _1 and a smaller pop- 
ulation at ~ 9 km s — 1 . Most of this revision was due to 
a more precise estimation of velocity uncertainties, but 
removing binaries did contribute to the reduction in a v . 

We did not obtain multi-epoch spectroscopy for 
VV124. Therefore, we were unable to remove the contri- 
bution of unresolved binaries. We suggest that the binary 
contribution to a v for VV124 is not nearly as large as for 
Bootes I or other ultra-faint dSphs. The warning s abou t 
the inflation of velocity dispe rsion bvlMinor et al. (l2010f) . 
IMcConnachie fc Cotel (|2010l) . and lKoposov et al.l ([20111 ) 
specifically concerned the ultra-faint dSphs with very low 
lumin osities (< 10 5 LgX V V124 does not fall into this 
class. Olszewski et al. (1996) showed that unresolved bi- 
naries inflate a v for the Ursa Minor and Draco dSphs by 
much less than 1 km s . Because the velocity dispersion 
of VV124 is closer to the dispersions of Ursa Minor and 
Draco than to ultra-faint dSphs, binarity likely does not 
dominate our measurement uncertainty. 

5.2. Rotation 

While IB 1 1 a! concluded that the H I gas in VV124 does 
not show any signature of rotation, the interpretation of 
the stellar wings as an edge-on disk of stars prompted us 
to search for stellar rotation. If the structure is a disk, 
it is nearly edge-on, and rotation in principle should be 
evident in measurements of stellar radial velocities. 

VV124 also contains H I gas that has a velocity gradi- 
ent of 5-10 km s _1 across the face of the galaxy and a full 
velocity range of ~ 45 km s . However , the g radient is 
not aligned with the disk-like structure. IBllal suggested 
that the gradient arises not from rotation, but asymmet- 
ric stellar winds acting preferentially on the southeast 
quadrant o f the g alaxy. 

Based on IB Hal s Figure 11 and our Figure [51 we define 
the "wings" of VV124 as \X\ > 4' (1.6 kpc) , where X 
is distance along the major axis. Our sample includes 
eight such member stars. Therefore, we have spectro- 
scopically confirmed the existence of th e phot ometrically 
feeble disk- like structure discovered by IBllal 

However, these stars show no evidence of rotation 
within the limitations of our measurements. Figure [8J 
shows stellar radial velocities versus distance from the 
minor axis. We quantified the amount of rotation as 
Awew, t ne difference between the average v r in the east- 
ern and western wings. Our sample contains six mem- 
bers in the eastern wing and two members in the western 
wing. The velocity measurements of these stars yielded 
Awew = 3.5 km s . 
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Table 3 

Properties of VV124 



Property 



Symbol 



Value 



Photometry 



1.36 ±0.03 Mpc 

8.2±{]l X 10 6 L e 

9A±ll x 10 6 M 

8.7 x 10 5 M & 

41.3" = 260 pc 



Distance 51 
Luminosity 15 
Stellar mass c 
Gas mass b 
Half-light radius b 



D 
W 

M* 



M, 



gas 
Re 



Dynamics 



Mean radial velocity 

Linc-of-sight velocity dispersion 

Mass within half-light radius d 

Mass-to-light ratio within half-light radius d 

Dynamical-to-stellar mass ratio within half-light radius°' d 

Total mass 

Total mass-to-light ratio 



(v r ) 

<J V 

M 1/2 

(M/L v ) 1/2 

(M dyn /M») 1/2 

Mtot 
(M/L v ) tot 



-29.1 ±1.3 km s" 1 

9.4 ±1.0 km s- 1 

(2.12 ±0.22) x 10 7 M® 

5.2 ±1.1 Mq/Lq 

4.5 ±1.9 

(1.95 ± 0.40) x 10 7 M Q 

4.8 ±1.3 Mq/Lq 



Metallicity 



Mean metallicity' 

Standard deviation 

Median metallicity 

Median absolute deviation 

Interquartile range 

Skewness 

Kurtosis 

Yield (Simple Model) 

Yield (Pre-Enriched Model) 

Initial metallicity (Pre-Enriched Model) 

Yield (Extra Gas Model) 

Extra Gas Parameter (Extra Gas Model) 



([Fc/H]) 
<r([Fe/H]) 



med([Fc/H 
mad ([Fe/H 
IQR([Fc/H 
Skew([Fe/H ) 
Kurt ([Fe/H]) 
p (Simple) 
o (Pre-Enriched) 
[Fe/H] 
p (Extra Gas) 
M 



-1.14 ±0.06 

0.49 
-1.08 

0.26 

0.54 
-0.64 ±0.31 
-0.05 ±0.60 



., n -+0.016 7 

n 117 +0.020 
U ' 1 - L ' -0.017 



'0 

z 



< -2.66 

0.109±°;°" Z Q 

4.76^; 79 



References. — a: Jacobs ct al. (20lJ). b: IBllal 
° Although Bll a] assumed M* / Ly = 2, we assumed that M*/Ly = 1.10, a value more typical of transition- 
type dwarfs (Woo ct al. 2008). We used this stellar mass-to-light ratio to calculate these quantities. 
The value could range from 0.8 to 1.5, and we took these limits into account in the estimation of the 

uncertainties. 

d Using the formula M 1/2 = 4G~ 1 iJ e at; l|Wolf et al.1120101 ). 

° Using the formula Mtot = 167/ir c a-J (Illingworth 1976), which may not be appropriate for VV124. This 

formula often gives smaller Mtot than the formula for M]/ 2 . 

Although Kirbv ct al. (2011a) calculated {[Fe/H]) with inverse variance weighting, we present the un- 
weighted mean. The low SNRs of many of our spectra cause the error on [Fe/H] to increase strongly with 
decreasing metallicity. Therefore, the weighted mean would be biased toward high [Fe/H]. 



In order to assess the precision of this measurement, 
we constructed theoretical rotation curves for the out- 
ermost regions of the galaxy. The curves were flat in 
the wings (r > 1.5'), and their magnitudes ranged from 
to 20 km s _1 , spaced at 1 km s _1 . We sampled the 
theoretical rotation curves with eight fake stars at the 
same distances along the major axis as the stars in our 
sample. Each fake star had a radial velocity perturbed 
from the rotation curve, and the amount of perturbation 
was sampled from a Gaussian random probability distri- 
bution with a standard deviation equal to the estimated 
uncertainty in v r for the corresponding observed star. 
Under the assumption that any rotation would be diluted 
by incoherent velocity dispersion, we further perturbed 
each star's velocity according to a Gaussian random dis- 
tribution with a standard deviation of a v , the previously 
measured line-of-sight velocity dispersion of VV124. We 
computed Awew for 1000 such Monte Carlo trials for 
each of the rotation curves. The standard deviation of 



the probability distributions is an estimate of the uncer- 
tainty on the observed value of Ai>ew = 3.5 km s _1 . This 
uncertainty, 8.6 km s , is insensitive to the magnitude 
of the rotation. 

Our measurements are consistent with the absence of 
rotation because Ai>ew is consistent with zero. How- 
ever, our data set is limited in its ability to detect rota- 
tion. Under the assumption that the velocity dispersion 
is constant with radius and that this velocity dispersion 
dilutes any rotation signal in the wings, the minimum 
rotation signal we could have detected is 8.6 km s _1 . 

The stars in the southeast of the galaxy do not have 
any more negative v r than stars elsewhere. Therefore, 
our observations do not show a stellar counterpart to the 
asymmetric gas velocity. Radial velocity measurements 
of bluer, younger stars might show such structure, but 
we have not performed such observations. 

Our observations are not conclusive evidence against 
rotation or a stellar counterpart to the H I velocity struc- 
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Figure 8. (a) Heliocentric radial velocity as a function of dis- 
placement from the minor axis (distance along the major axis). 
Filled points are spectroscopic members, and hollow diamonds and 
crosses are non-members. The horizontal dashed line indicates 
(v r ). Only stars passing all membership cuts except velocity are 
shown, (b) Spatial locations of members and stellar non-members 
fainter than the TR GB. The gray lines are contours of equal sur- 
face density (BllaD. From the inside to the outside, the contours 
represent 3<r, 4<r, 5<r, 6<r, 7<r, 8<r, and 9cr differences from the central 
surface density. 

ture. Furthermore, we may have misclassified some fore- 
ground stars as members. Such misclassifications may 
alter a rotation signal. We considered the possibility that 
some true members may have been misclassified as fore- 
ground dwarfs based on the strength of their Na I 8190 
doublets (hollow diamonds in Figure [3]) , but the radial 
velocities of those stars do not suggest rotation or veloc- 
ity structure any more than the spectroscopic members. 
We cannot rule out the possibility that a larger sample 
with more precise measurements of v r will detect velocity 
structure in the stars. 

5.3. Mass 

The mass of a spheroidal stellar system can be deter- 
mined from its line-of-sight velocity dispersion under the 
assumption of a mass profile and velocity anisotropy. In 
the case of VV124, only stars and dark matter signifi- 
cantly affect cr v . Gas is irrelevant because t he gas mass 
is just a small fraction of the stellar mass. Illingworth 
(|1976f ) presented a formula to measure the mass of glob- 
ular clusters: M tot = 167^r c cr^, where r c is the iKingl 
(1962) core radius and /Xjis a dimensionless number that 
quantifies the mass concentration toward the system's 
center. The core radius may be approximated in rela- 
tion to R e , t he half - light or effective radius, as follows: 
r c = 0.64i? e . iMateol (| 19981 ) stated that \x = 8 is appropri- 
ate for low-concentration King profiles and therefore ap- 



7 IKi ng (1966) and Illingworth 119761) called this parameter fi. 
IMateol 11998) and lSimon fc G cha (2003) called it /3. Because is 

typically used to represent velocity dispersion anisotropy, we have 
chosen to call the concentration parameter fj,. 



JTopriate for most dwarf galaxies. However. Illlingworthl 's 
19761) formula assumes that mass follows light and that 
the velocity dispersion anisotropy j3 = 0. These two as- 
sumptions are appropriate for globul ar clusters but not 
for dark matter-dominated galaxies. IWolf et ail (2010) 
presented an alternative formula for the mass enclosed 
within the half-light radius: Miy 2 = 4G _1 i? e cr^, where 
R e is the half-light or effective radius. It is at this ra- 
dius that /3 has the least influence on the estimate of 
the enclosed mass. Inside and beyond R e , the enclosed 
mass becomes increasingly uncertain because /3 is un- 
known when only line-of-sight ve locities are k nown . Ta- 
ble GH gives both M tot fr om the llllingworthl (|1976[ ) for- 
mula and Mi/2 from the IWolf et all (|2010D formula. As 
is typical for dwarf galaxies, Mi/ 2 actually exceeds M tot 
because the llllingworthl (|1976l ) formula is inappropriate 
for dwarf galaxies. Note that both of these formulae 
assum e sphericity. The ellipticity of VV124 (e = 0.44, 
IBllal ) diminishes the accuracy of these mass estimators 
for VV124. 

The mass-to-light ratio within the half-light radius is 
(M/L v ) 1/2 = 5.2 ± 1.1 M q /Lq. This value is too 
large for a mass dominated by stars alone. Accord - 
ing to up-to-date stellar models (jPercival et al.l [2009) . 
old (12 Gyr) and metal-poor ([Fe/H] = —1, which is 
about the mean metallicity for VV124, as derived in Sec- 
tion [B]) stellar populations without dark matter exhibit 
M/Ly < 2.5 Mq/Lq. Due to the contribution of the 
young population, the expected value for VV124 should 
be lower. The ratio of the dynamical mass to the stel- 
lar mass within the half-light radius is {M < x ya /M if ) 1 /2 = 
4.5 ± 1.9. (This value depends inversely on the assumed 
stellar mass-to-light ratio, M^jLy. See footnote c in Ta- 
ble [3]) Therefore, we conclude that the mass of VV124 
is domi nated by dark matter. 

IBllal roughly estimated the mass-to-light ratio 
(M/L v w 8) of VV124 from the H I velocity distribution. 
They assumed that the gravitational potential dominates 
the gas motion and that the H I distribution is spherical 
and isotropic. Their value is 1.5 times our measurement 
of (M/Lv)i/2- However, theirs was a rough estimate. 
The difference can be explained by their slightly larger 
measurement of 11 km s _1 for the velocity dispersion, 
compared to our measurement of 9.4 km s . Our mea- 
surement is more accurate because it is based on individ- 
ual stars. The gas velocity dispersion could be sensitiv e 
to stellar feedback, especially in light of Bella zzini et alj 's 
interpretation of stellar winds as the cause of the velocity 
gradient seen in the gas. 

6. METALLICITY DISTRIBUTION 

We measured iro n abundance s from Fe I l i nes in the 
DEIMOS spectra. IKirbv et all (|200l 120091 l2010l) de- 
scribed the procedure that we used for these measure- 
ments. We cor re cted t he photometry for extinction with 
ISchlegel et all 's (|1998l) dust maps. The average redden- 
ing was low (E(B — V) ~ 0.01). After we converted 
the extinction-corrected apparent magnitudes to abso- 
lute magnitudes as suming (m — M) = = 25.67 ± 0.04 
( 
( 



Jacobs et al. 



Girardi et al 



20111) we used 12 Gyr Padova isochrones 
2002) to estimate the photometric effec- 
tive temperature and surface gravity of each star. A 
large grid of synthetic spectra was searched for the best- 
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Figure 9. Metallicity measurement differences of the same stars 
observed on two different siitmasks. The differences are divided 
by the quadrature sum of the measurement uncertainties, which 
include noise and systematic error terms. The solid curve is a 
Gaussian with unit standard deviation. The systematic error term, 
°sys,[Fe/H]! was adjusted to ensure that the standard deviation is 
unity. Compare this figure to Figure [4] 
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Figure 10. Top: Absolute metallicity differences of repeat mea- 
surements of the same stars on two different siitmasks versus SNR. 
These measurements were necessarily based on spectra with about 
half the exposure time of the final, stacked spectra. Bottom: Esti- 
mated error on [Fe/H] as a function of SNR for all member stars. 
The total error is the quadrature sum of the random and systematic 
error terms. Compare this figure to Figure [5] 



fitting synthetic spectrum. In the search, temperature 
and [Fe/H] were varied. The temperature was assigned 
a probability distribution based on the uncertainty in 



the star's position in the CMD. This probability distri- 
bution was included in the x 2 minimization to find the 
best-fitting synthetic spectrum. In this way, the spec- 
troscopic temperature was not allowed to stray far from 
th e photometr i c tem perature. (See Sections 4.5 and 4.7 
of IKirbv"et al. 2009 for more details.) The photometry 
and the spectroscopy shared roughly equal weight in de- 
termining the final temperature. The photometry alone 
determined the surface gravity. 

In previous works based on this technique, the un- 
certainty in [Fc/H] (cr covar ) was calculated in part from 
the covariance matrix in the Levenberg-Marquardt op- 
timization. Various diagnostics indicated that this ran- 
dom uncertainty term underestimated the full error in 
the measurement. Therefore, a systematic error term, 
"sysJFe/H] = 0.113, was added in quadrature to the ran- 
dom uncertainty. That term was calculated from spectra 
with a typical exposure time of one-hour. The exposure 
time for most spectra in this work exceeds seven hours. 
Therefore, errors in sky subtraction play a much larger 
role than in our previous works. 

Because the spectral noise behavior assumes a differ- 
ent character than in the spectra from which the system- 
atic error term was originally calculated, we re-evaluated 
the magnitude of cr S ys,[Fe/H] using only the VV124 data 
set. In analogy to our calculation of radial velocity un- 
certainties, we examined differences in [Fe/H] between 
measurements of the same stars on the two separate slit- 
masks. We determined cr S ys,[Fc/H] — 0.191 from the fol- 
lowing equation. 



\ 



1 N 

-Y 



([Fe/H],- [Fe/H] 2 



x ^covar,! + a 



2 
covar,2 



2cr sys,[Fc/H] 



= 1 



(2) 



In other words, a. 



sys,[Fc/H] 



0.191 is the value required to 



force la agreement between the repeat measurements of 
[Fe/H]. Figure IH1 shows that this agreement was achieved. 

Figure [TO] shows the relationship between the error or 
uncertainty on metallicity versus spectral SNR. The top 
panel shows the differences in [Fe/H] measured from the 
two different siitmasks. These measurements are dif- 
ferent by up to 0.8 dex because most of them (those 
observed on both siitmasks) are based on partial expo- 
sures with low SNRs. The bottom panel shows the es- 
timated uncertainty on metallicity versus SNR. These 
estimates are based on the full exposures (coadded spec- 
tra from both siitmasks). Therefore, the bottom panel 
shows larger SNRs than the top panel. The trend of 
error with SNR is weaker than for radial velocities (Fig- 
ure [5]) because cmc for radial velocities was calculated 
with Monte Carlo resampling of the observed spectra. 
Monte Carlo resampling was computationally prohibitive 
for estimating uncertainties on [Fe/H]. As a result, we re- 
lied in part on a covar instead of iJmc,[Fc/h] for estimating 
the uncertainty on [Fe/H] . This quantity comes from the 
covariance matrix, which tends to underestimate uncer- 
tainty compared to Monte Carlo techniques. Therefore, 
u sys rpe/H] accounts for the majority of the uncertainty 
in most stars, and the resulting error distribution shows 
only a weak trend with SNR. 

In the following discussion, we have assumed that 
the solar iron abundance is 12 + logn(Fe)/ra(H) = 7.52 
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Figure 11. Metallicities derived from the CaT versus metallicities 
derived from spectral synthesis. 



(|Sneden et al.f J992). Also, we discarded metallicity mea- 
surements with uncertainties larger than 0.5 dex. 

Another popular method of measuring metallicities of 
red giants from low- to medium-resolution spectra is the 
empirical relation between [Fe/H] and the EW of the 
CaT. In order to compare CaT metallicities to our syn- 
thetic spectral metallicities, we measured the EWs of the 
CaT lines at 8542 A and 8662 A. We fitted Lorentzian 
functions to each of the CaT lines. We estimated un- 
certainties on EWs by Monte Carlo resampling. We 
added Gaussian random noise to each spectrum in pro- 
portion to its error spectrum. The amount of noise 
added to each pixel was sampled from a Gaussian ran- 
dom distribution with a standard deviation equal to that 
pixel's estimated Poisson noise. We took the error on 
EW as the stand ard deviation of 10 such realizations. 
Then, we applied IStarkenburg et all 's (|2010l) calibration 
between CaT EW, absolute V magnitude, and [Fe/H]. 
This calibration was based on a combination of Gaus- 
sian fits and numerical integration in order to include 
the pressure-broadened wings of the CaT. As a rough 
substitute, we used the Lorentzian, which is a good ap- 
proximation to both the cores and the wings of the CaT 
lines. We transformed extinction-correct ed go and r p 
magnitudes to Cousins Vq magnitudes using lJordi et all s 
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Figure 12. Metallicity distribution for member stars. Error bars 
are calculated with Poisson statistics. Best-fitting, analytic models 
of chemical evolution are shown as colored curves. The red and blue 
lines overlap almost exactly. 



(2006) equations. We computed errors on [Fe/H] CaT by 
propagating the errors on EW and V. We did not in- 
clude the uncertainty in the metallicity calibration be- 
cause it is significantly less than the substantial errors 
in EW. Figure [TT] shows the comparison between [Fe/H] 
measured from the CaT and from spectral synthesis for 
the stars where it was possible to measure Lorentzian fits 
confidently. Although the measurement uncertainties are 
large, the agreement is good. This comparison reinforces 
credibility in our metallicity measurements based on Fe I 
lines. For the rest of this discussion, we used synthetic 
metallicities instead of the CaT metallicities because the 
synthetic metallicities are more precise for two reasons. 
First, the synthetic metallicities do not require a conver- 
sion between photometric systems. Second, Fe I lines in 
the DEIMOS spectrum of a typical moderately metal- 
poor red giant have more signal (sum of EWs) than the 
combined CaT. 

Figure [T^] shows the metallicity distribution of mem- 
ber stars. The metallicities range from [Fe/H] = —2.3 
to +0.0. The shape of the distribution (a peak with a 
longer metal-poor tail than a metal-rich on e) is typical 
of dwa rf galaxies of VV124's stellar mass (|Kirbv et al.l 
I2011al) . The bottom of Table [3] lists some additional 
shape characteristics of the metallicity distribution. 

We measured a mean metallicity of ([Fc/H]) = 
— 1.14 ± 0.06. This number is larger than prc- 
yious photometric estimates: [Fe/H] = —1.37 
(IKopvIov et al.l [2001: iTikhonov et all [20101 ) . [Fe/H] = 
-1.79 II Jacobs et al.l [2011J), [M/H] = -1.5 (in the 
outskirts. IB Hal ), and [M/H] = —1.4 (in the region 
0.5' < R f < 1 .0', wh ere R € is the elliptical radius, 
IBellazzini et al.l 12011b ). The metallicities measured 



bv I Jacobs et al.l (|20 111) and IBellazzini et al.l (|2011aHb1 ) 
are con sistent because I Jacobs et al.l mea sured [M/H] 
whereas iBellazzini et al.l measured [Fe/H] . Uacobs et al.l 



12 



Kirby, Cohen, & Bellazzini 



(|2011l) and lBellazzini et al.l (|2011allbh also noted that the 
g reat majority of old RGB stars have [Fe/H] < — 1. 

iKirbv et al.l ()2011al ) fitted three analytic chemical evo- 
lution models to the metallicity distributions of eight 
dSph satellites of the Milky Way. The Simple Mode l 
is also known as a Closed Box (Tal bot fc Arnettlll971h . 
Its single parameter is the effective yield, p. Low effec- 
tive yield indicates that the galaxy lost metals. Hence, 
we prefer the name "Simpl e Model" to "Closed Box." 
The Pre-Enriched Model (|Pagell 119971 ) is the Simple 
Model with the additional parameter [Fc/H] , which is 
the initial metallic ity of the gas. The Extra Gas Model 
(|Lvnden-BeUlll975J) allows for gas to fall into the galaxy 
during star formation. The infall rate is specifically cho- 
sen to allow for an analytic solution to the metallicity 
distribution. The parameter M quantifies the amount 
of infalling gas. The Extra Gas Model reduces to the 
Simple Model in the limit M = 1. 

We fitted the same models with the same maximum 
likelihood estimation procedure as IKirbv et al.l ()2011aD . 
The colored curves in Figure [12] show the results. The 
curves for the Simple and Pre-Enriched Models nearly 
overlap because the Pre-Enriched Model reduces to the 
Simple Model in the limit that [Fe/H] = — oo. The max- 
imum likelihood estimation could not constrain a lower 
bound to [Fe/H] . Hence, the metallicity distribution 
is consistent with initially metal-free gas. The metallic- 
ity distribution is more peaked than the Simple Model. 
The Extra Gas Model roughly accounts for this peak. 
In fact, the Extra Gas Model is 5.9 times more likely 
to describe the metallicity distribution than the Simple 
Model. Metallicity distributions consistent with gas in- 
flow are typical for d Sphs with the stellar mass of VV124 
(|Kirbv et al.ll2011aD . 

The low metallicity and highly sub-solar effective yields 
(p ~ 0.1 Ztn) require th at VV124 lost some of its metals 
(see IKirbv et al.l 1201 lbl) . possibly in the form of super- 
nova ^whids7[BTla|also noted the role of gas flows into 
and out of VV124. The H I gas distribution is asymmet- 
ric and suggestive of gas flowing to or from the southeast 
of the galaxy's star-forming center. The densest con- 
centration of gas is also offset to the southeast, further 
supporting the bulk motion of gas along that direction. 
While the present gas motions did not affect the past 
star formation that made the bulk of the RGB, the radio 
observations do support the role of gas flows in shaping 
the stellar population of VV124. 

Radial metallicity gradients are typical of many 
galaxies, a nd many dwarf g a laxies are no excep 
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Koch et al. 2006; Spolaor et al 



Palla dia et alj 120061: 
[Koleva et al.l 

2009a, b; Kirb v et alll2011al ). If gradients exist, they are 
always negative in the sense that metallicity decreases 
with increasing radius. Photometry strongly indicated 
that the stars near the center of VV124 are younger 
and more metal-rich (Tikho nov et al|[2010t|jacobs et al.l 
120111 : Bellazzini et al. 2011b]). We confirmed a gradient 
of d[Fe/H]/dr = -0.26 ± 0.04 dex kpc" 1 . Figure Q3 
shows the radial distribution of [Fe/H] along with an 
illustration of the gradient. Radial gradients can arise 
by several mechanisms. The most plausible explanation 
is that metallicity increased with time as supernovae en- 
riched the interstellar medium with iron. Star formation 
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Figure 13. Metallicity as a function of elliptical distance from 
the center of VV124. The elliptical distance is the semi-major axis 
of the ellipse (e = 0.44) on which the star lies. The red line is a 
least-squares fit. 

also became more centrally concentrated because gas loss 
(due to supernova winds, for example) preferentially ex- 
pelled gas from the outer reaches of the galaxy. There- 
fore, the radial gradient arose because successive gener- 
ations of stars became simultaneously more metal-rich 
and more centrally concentrated. 

7. DISCUSSION AND CONCLUSIONS 

Since its recent relocation to the edge of the Local 
Group, the dwarf galaxy VV124 has received consider- 
able scrutiny. However, measurements of its radial ve- 
locity have been inconsistent. The galaxy also seemed 
to possess a disk of stars. We observed VV124 with 
Keck/DEIMOS, and we determined that its radial ve- 
locity agrees with H I measurements. The stellar veloc- 
ities we measured show no evidence of a disk rotating 
any faster than 8.6 km s _1 , but larger samples will be 
necessary to place stricter limits on rotation. 

We gave an estimate of the dynamical mass of VV124 
and determined that the mass of its stars alone does not 
account for the velocity dispersion. Therefore, VV124 
contains dark matter. Its mass-to-light ratio agrees with 
dSphs of similar stellar mass. The metallicity distribu- 
tion is also similar to dSphs of similar stellar mass. The 
average metallicity, metallicity dispersion, and shape of 
the distribution are nearly indistinguishable from dSphs 
like Leo I or Fornax. As is typical for large dSphs, the 
metallicity distribution fits a chemical evolution model 
with infalling gas better than a Closed Box Model. 

The parameters we derived from new spectroscopy sug- 
gest that VV124 is a dlrr far along on its transformation 
into a dSph. Figure 1141 shows that VV124 lies along the 
same trends defined by other dSphs and dwarf elliptical 
galaxies in the Local Group. Thus, the galaxy is an ex- 
cellent test case for models that posit that dlrrs evolve 
into dSphs by the removal of gas and randomization of 
stellar velocities. This process may happen in isolation 
or in the vicinity of a much larger galaxy, such as the 
Milky Way. 

The transformation of a dlrr into a dSph might pro- 
ceed without the benefit of a dense environment. For 
example, VV124 in its present form could have resulted 
from a merger of two dwarf galaxies ~ 1 G yr ago 
(|D'Onghia et all 120091: [Kazantzidis et alj|2011b|) . The 
merger would have disrupted a disk and possibly trig- 
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Figure 14. Velocity dispersion, mass within the half-light radius, mass-to-light ratio within the half-light radius, dynamical-to-stellar 
mass ratio within the half-light radius, and average metallicity versus stellar mass (left) and luminosity (right) for dSph and dwarf elliptical 
galaxies in the Local Group. VV124 is represented by a red star. Dynamical quantities (<r v , Afi/2, (M/Lv)x/2i an d (A^dyn/-^*)l/2) were 
taken from [Wolf ct al. (201CJ) and references therein for all galaxies except Boo tes I. For Bootes I, we adopted cr v from lKoposov et alJ 
(2011), and we adjusted its mass accordingly. The stellar masses were taken from Woo et al. (2008) for the larger dSphs and Martin et al. 
(120081, using the va lues derived with t he [Kr oupa et al. 1993 initial mas s function) for the ultra-faint dSphs. The metallicities came from 
Kirby ct al. (2011a) except for Carina (Hclmi ct al. 2006) and Bootes I (Martin ct al. 2007). 



gered the starburst that Uacobs et all (J2011D observed. 
Recent simulations showed that ultraviolet background 
radiation and stellar feedback are sufficient to exp el all 
of the gas from a dwarf galaxy (iSa wala et alj |2010). The 
challenge in these models is the formation of disk-less 
dSphs. It is possible that VV124 could have formed as a 
dispersion-supported galaxy, but this scenario runs con- 
trary to the view that dwarf galaxies formed as disks from 



which super nova feedback removed gas with low angular 
momentum (|Governato et al.l I2010T ). Another challenge 
is to explain the low gas fraction of VV124 relative to 
other isolated dwarf galaxies. In a sample o f 39 d warf 
galaxies with ste llar masses fr o m I Woo et al.l (2008) and 
gas masses from IGrebel et al.l ()2003D . the only galaxies 
with lower gas fractions than VV124 (/ gas = 0.09) are 
satellites of the Milky Way or M31, not isolated dwarfs 
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like VV124. 

The tidal stirring model (jMaver et all l2001allbl ) ex- 
plains the transformation from dlrrs to dSphs in the con- 
text of galaxy environment. This model posits that disky 
dlrrs lose their gas through ram pressure stripping in a 
group environment. The stars then lose their rotational 
structure by tidal interactions. Numerous simulations 
address the transf ormation of dlrrs to dSphs by tidal stir- 
ring (e.g., Klimcntows ki et ail 12009a] ; lLokas et al.l 120101 
l2011tlKazantzidis et alJl2011al ). The stated goal of most 
of these simulations is to construct a spheroid from a 
disk. As a result, few studies address the kind of mor- 
phology peculiar to VV124 (a spheroid with a very faint 
disk- like feature). Nonetheless, the intermediate stages 
of some tidal stirring simulations resemble a spheroid 
with a disk. The simulated tidally stirred dwarf galaxy 
in Figure 12 of lMayer et al.l (|2001b[ ) is a decent represen- 
tation of VV124, but the "disk" is actually tidal tails. 
Our radial velocity observations of stars in the wings 
of VV124 do not support their interpretation as a disk. 
Hence, they may be tidal tails. However, tidal tails near 
apocen ter are expected to be elo ngated toward the per- 
turber (Klimcntowski et al. 2009b) . The disk- or tail- like 
feature observed in VV124 must be nearly perpendicu- 
lar to the Milky Way, and the great circle that connects 
VV124 and M31 intersects the major axis of VV124 at 
an angle of 54°. Therefore, the elongated photometric 
features are unlikely to be tidal tails resulting from an 
interaction with either of the large spiral galaxies of the 
Local Group. 

Furthermore, VV124 is in a remote corner of the 
Local Group, highly is olated from other galaxies. 
iKarachentsev et all (|2009() concluded that VV124 resides 
outside the zero- velocity surface of the Local Group. Be- 
cause its velocity relative to the Local Group's barycen- 
ter is positive, it will ne ver enter the group in the fu- 
ture. Additionally, iBllal calculated that the intragroup 
medium around VV124 is probably not dense enough to 
remove gas by ram pressure stripping. As a result, tidal 
stirring would require that VV124 passed near another 
galaxy and subsequently traveled ~ 1 Mpc to its current 
position. Traveling such a large distance after a tidal 
interaction is not out of the question. A dwarf galaxy 
that passed within half of the virial radius of a mas- 
sive halo could have trav eled up to several virial radii 
beyond that massive halo ([Lin et al.ll2003t iLudlow et al.l 
I2009L iTevssier et all 1201 21 ). As many as 10% of all dark 
matter halos are actually ejected subhalos that can be 
found as far away as fo ur virial radii from their former 
hosts (|Wang et alj |2009). VV124 qualifies as a candidate 
subhalo ejected from either the Milky Way or M31. Its 
low velocity relative to the Local Group would require 
VV124 to be near its apocenter. Such a location would 
not be surprising given its large distance from the nearest 
massive halo. It also makes sense that we would observe 
an ejected subhalo near its apocenter because subhalos 
linger near their apocenters. As a potentially ejected 
satellite galaxy, VV124 could be an example of tidal stir- 
ring. In this scenario, a brief, high-speed passage near 
the Milky Way or M31 removed most of its gas, and the 
violent change in tidal force near pericenter disrupted 
its stellar disk. A single passage is sufficient for ram 
pressure stripping to remove most o r all of the gas in a 
galaxy of VV124's circular velocity ([Maver et al.ll2006h . 



but most simulations of tidal stirring require multiple 
pericentric passages to t ransform a disk into a spheroid 
i|Kazantzidis et al.ll2011af) . 

The tidal stirring explanation requires special circum- 
stances. First, VV124 had to have been ejected from the 
Milky Way or M31, possibly in a multi-body encounter. 
Second, the stellar disk would have to be disrupted in just 
one pericentric passage, or VV124 would have had to ac- 
quire a large amount of orbital energy in a multi-body 
encounter on a subsequent pericentric passage. Third, 
the tidal tails — if indeed the photometric wings are tidal 
tails — would have to somehow point in a direction away 
from the former host galaxy. 

While our observations support the interpretation of 
VV124 as a galaxy of a type between dlrr and dSph, 
they do not explain how it came to be that way. We 
look forward to future simulations that will manufacture 
a digital replica of VV124: an isolated, mostly spheroidal 
dwarf galaxy that somehow lost most of its gas but re- 
tained wispy, minimally rotating wings. 
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Table 2 

Target List 



ID RA (J2000) Dec (J2000) 



Masks a SNR (A" 1 ) v r (km s" 1 ) EW(Na I 8190) (A) 



[Fe/H] 



Member' 



Reason 



'< 

O 

c 

a 
p 

a 



13585 


09 15 09.6 


+52 49 52.0 


24.351 ± 0.030 


23.295 ± 0.043 


1 


1.1 










N 


G 


16765 


09 15 10.8 


+52 48 34.1 


25.600 + 0.038 


23.667 ± 0.046 


1 


2.9 










N 


Bad 


11688 


09 15 11.8 


+52 50 12.4 


25.268 + 0.029 


23.365 ± 0.025 


1 


15.5 










N 


G 


16355 


09 15 11.8 


+52 48 50.4 


24.183 + 0.011 


23.125 + 0.016 


1 


3.5 










N 


Bad 


16895 


09 15 13.4 


+52 48 29.4 


24.198 + 0.021 


23.337 ± 0.032 


1 


2.9 










N 


Bad 


3705 


09 15 15.4 


+52 51 59.3 


24.651 + 0.015 


23.217 + 0.018 


1 


25.6 


+31.5 + 2.7 


3.21 + 0.24 






N 


v r Na 


16598 


09 15 16.9 


+52 48 40.7 


23.009 ± 0.006 


21.553 + 0.007 


1 


50.9 


-68.1 + 2.4 


2.22 + 0.06 






N 


v r Na Bright 


16798 


09 15 17.5 


+52 48 32.7 


22.361 ± 0.006 


20.841 ± 0.007 


1 


68.8 


-25.5 + 2.3 


2.45 + 0.04 






N 


Na Bright 


7035 


09 15 17.9 


+52 50 58.8 


22.720 ± 0.010 


22.720 + 0.013 


2 


16.5 










N 


G 


7917 


09 15 19.8 


+52 50 49.5 


24.287 + 0.012 


23.183 + 0.017 


2 


13.0 


-40.1 + 4.8 




-1.02 + 0.21 


Y 





References. — Identifications, coordinates, and photometry from Blla. 

Note. — (This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 
a Number of DEIMOS masks on which the object was observed. 

b Reasons for non-membership. v r : Inappropriate radial velocity. Na: Spectrum shows strong Na I A8190 doublet. G: Spectrum shows emission lines or redshifted Ca H and K 
indicating that the object is a galaxy. Bright: Target is brighter than the TRGB. Bad: Spectral quality was insufficient for radial velocity measurement. 



lines, 



